Meta-nematic transitions in a bilayer system: Application to the bilayer ruthenate 
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It was suggested that the two consecutive metamagnetic transitions and the large residual resis- 
tivity discovered in SrsRu207 can be understood via the nematic order and its domains in a single 
layer system. However, a recently reported anisotropy between two longitudinal resistivities induced 
by tilting the magnetic field away from the c axis cannot be explained within the single layer nematic 
picture. To fill the gap in our understanding within the nematic order scenario, we investigate the 
effects of bilayer coupling and in-plane magnetic field on the electronic nematic phases in a bilayer 
system. We propose that the in-plane magnetic field in the bilayer system modifies the energetics 
of the domain formation, since it breaks the degeneracy of two different nematic orientations. Thus 
the system reveals a pure nematic phase with a resistivity anisotropy in the presence of an in-plane 
magnetic field. In addition to the nematic phase, the bilayer coupling opens a different route to a 
hidden nematic phase that preserves the x-y symmetry of the Fermi surfaces. 
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I. INTRODUCTION 

Electronic liquid crystal phases have been widely dis- 
cussed in the context of doped Mott insulators^ and 
have attracted much attention with the discovery of 
anisotropic quantum Hall phases in GaAs heterostruc- 
tures in large magnetic fields 4^ In particular, the 
anisotropy in the longitudinal magnetoresistivities in the 
bilayer ruthenate Sr3Ru2 07 strongly suggests the exis- 
tence of an anisotropic metallic phase dubbed an elec- 
tronic nematic phased 

In experiments on ultrapure Sr3Ru207 in magnetic 
fields along the c axis, an unusual phase characterized 
by a pronounced residual resistivity emerges in the vicin- 
ity of a putative quantum critical point. Furthermore, 
this phase is bounded by two consecutive meta-magnetic 
transitions ASaSftiSftH Upon tilting the field slightly to- 
wards one of the in-plane crystal axes, a magnetoresis- 
tive anisotropy appears, where the pronounced resistive 
anomaly parallel to the in-plane field direction remains 
unchanged, but disappears in the perpendicular direc- 
tion. 

It was proposed that the two consecutive meta- 
magnetic transitions occur due to the formation of ne- 
matic order As a consequence of the nematic order, 
one expects to find an anisotropy in the longitudinal re- 
sistivities in the pure nematic phase due to the Fermi 
surface distortion, but this has not been observed in 
Sr3Ru207 in magnetic fields along the c axis. Instead, 
the pronounced resistivity shows up in the proposed ne- 
matic region bounded by two metamagnetic transitions. 
The large resistivity was explained by scattering due to 
domains of two degenerate nematic orientations, 13 How- 
ever, when the magnetic field is tilted away from the c 
axis, the transport anisotropy is discovered. This cannot 
be understood within the previous nematic order pro- 
posal, since the Zeeman coupling and the energetics of 
the domains are independent of the magnetic field direc- 
tion in the single layer system. 



In this paper, we attempt to understand the recently 
reported magnetoresistive anisotropy in the presence of 
an in-plane magnetic field. We show that the magnetore- 
sistive anisotropy in the presence of an in-plane magnetic 
field can be understood within the nematic order pic- 
ture when one takes into account the bilayer coupling. 
Note that the in-plane magnetic field in the bilayer lat- 
tice causes a relative momentum mismatch between the 
layers. When the in-plane field is along one of the crystal 
axes, it breaks the degeneracy of two different nematic 
orientations. Consequently, domains with different ne- 
matic orientations are no longer energetically favorable, 
and the system exhibits a pure nematic phase with a 
Fermi surface elongation. Thus, the anisotropy in trans- 
port is recovered in the presence of an in-plane magnetic 
field. 

The paper is organized as follows. We introduce the 
bilayer model in Sec. II. In Sec. Ill, we identify the dis- 
tinctly different nematic phases and present the phase 
diagram as a function of the bilayer coupling and chem- 
ical potential. We also discuss a hidden nematic phase 
which is absent in the single layer system. An in-plane 
magnetic field is incorporated in Sec. IV, where we study 
the phase diagram under the in-plane magnetic field and 
the signatures of meta-nematic transitions in the longi- 
tudinal conductivity and the magnetic susceptibility. We 
close with a discussion and a summary of our findings in 
Sec. V. 



II. BILAYER MODEL 

In the electronic nematic phase, electron momenta pre- 
fer to be aligned along a certain direction, typically along 
one of the crystal axes, thus breaking a point-group sym- 
metry of the underlying lattice. A number of models 
have been employed to study the formation of electronic 
nematic order. The approach adopted in the present 
work is based on the idea of a broken symmetry state 
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of an isotropic liquid ^^^^i^^^ 
Here, the formation of the nematic phase is due to a 
spontaneous Fermi surface distortion often referred to as 
a Pomeranchuk instability. 28 It was found, however, that 
the divergence of the nematic susceptibility, which defines 
the Pomeranchuk instability, is preempted by a first or- 
der transition, and that the formation of an electronic 
nematic phase on a lattice is intimately connected to the 
van Hove singularities in the density of states (DOS) i 14 i 15 
The effective nematic interaction successfully describes 
several novel phenomena observed in S^R^Oy*^^ We 
generalize the quadrupole density interaction for a bilayer 
square lattice as follows, 
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Here, A = ±1 denotes the layer index, a = ±1 the spin 
degree of freedom, and c\ c the electronic creation and 
annihilation operators. The functions ^(q) an d ^(q) 
denote intraplane and interplane quadrupolar density in- 
teractions, respectively. The order parameter is defined 
through the tensors (Q^) in analogy to its counterpart 
in classical liquid crystal theory. In conjunction with the 
tight-binding model on a square lattice, the intralayer in- 
teraction term in Eq. (pQ) describes the first order transi- 
tion between isotropic and nematic states within a mean- 
field theory,^ 

In general, there are two distinct nematic orders in 
a single layer square lattice system. The preferred di- 
rection of the electron momenta can be aligned either 
parallel or diagonal to the crystal axes. Previous studies 
of monolayer systems have shown, however, that diag- 
onal order is generally suppressed for the model given 
by Eq. ([T]) j 29 i 3Q i 31 We, therefore, concentrate on the ne- 
matic phase parallel to the in-plane crystal axes. Assum- 
ing an attractive interlayer interaction potential F2 (q) = 
—F2 £ q ,o (i*2 > 0) and a generally attractive intralayer 
interaction ^(q) = — G2#q,o> the components of the 
parallel order parameter are defined by 
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where a positive (negative) value signals, that electron 
momenta are preferentially aligned along the y (x) axis. 
Note that the order parameter is defined layerwise (A = 
±1). Within the mean- field (MF) approximation, the 
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where t± is the interlayer hopping amplitude. 

To consider the effect of an in-plane magnetic field, 
we incorporate the field via Peierls substitution and Zee- 
man coupling. Assuming that B = B x x, the magnetic 
flux in the x-direction through a single plaquette is given 
by (j) x = B x a 2 , as the layer separation in Sr3Ru207 is 
of about the same order as the planar lattice parame- 
ter a. Normalized by the flux quantum O = hc/e, the 
flux causes the following relative momentum mismatch 
between the layers: 
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The electronic dispersions for each layer in the presence 
of an in-plane magnetic field are given by 
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where k — > k + p for the lower layer (A = — 1). In Eq. 
(j6]), we introduced the chemical potential fi and an ef- 
fective magnetic moment 7 (= gfiB-^j), while a = ±1 
stands for spin-up and spin-down. Taking into account 
the bilayer coupling t± , the hybridized energy bands read 
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where the order parameter components satisfy the set of 
mean- field equations (A = ±1) 
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The expectation values depend in a non-linear fashion on 
all . Both dispersions of Eq. ([7j) are mapped onto 
each other under a particle-hole transformation. 

In the following, we calculate the free energy using an 
adaptive integration scheme and solve the set of mean- 
field equations self-consistently. We first present the ef- 
fects of the bilayer couplings t± and G2 on nematic order 
in the absence of a magnetic field (cj) x = 0). Based on the 
zero-field results, we then deduce the effect of an in-plane 
magnetic field originating from the bilayer coupling t±. 
Finally, we include the Zeeman term and investigate the 
signatures of nematic order in magnetization and trans- 
port. 




FIG. 1: (Color online) Fermi surface shapes of the two ne- 
matic phases in the bilayer system. The solid lines represent 
the prehybridized and the dotted lines the hybridized Fermi 
surfaces. The left panel shows the Fermi surface configura- 
tion in the parallel nematic phase, where both prehybridized 
Fermi surfaces are elongated parallel to each other. The right 
panel shows the hidden nematic phase, where both prehy- 
bridized Fermi surfaces are elongated perpendicular to each 
other, while the hybridized Fermi surfaces recover the fourfold 
symmetry. 



III. NEMATIC PHASES IN A BILAYER 
SYSTEM 

Let us first identify possible nematic phases in the ab- 
sence of external fields (<p x = 0). The spin index is redun- 
dant and the state of the system is specified by a layer- 
dependent order parameter A^ A ^ (A = ±1). As none 
of the crystal directions is distinguished, both compo- 
nents of the order parameter have the same magnitude, 
|AW| = |A( -1 )|. This implies that, besides the isotropic 
phase, where A^) = A^ -1 ) = 0, only two^ distinct ne- 
matic phases can occur as shown in Fig. 1: a paral- 
lel nematic phase, where A^) = A^ -1 ) and a "hidden" 
nematic phase, where A^ 1 ^ = — A^ _1 \ In the hidden 
nematic phase, the prehybridized Fermi surface of one 
layer is elongated along the x (y) direction, while 

the prehybridized Fermi surface of the other layer 
is stretched along the y (x) axis. However, taking into 
account the bilayer coupling, each of the hybridized en- 
ergy bands shown as the thick lines in Fig. [T] preserves 
the x-y symmetry, but breaks the relative symmetry be- 
tween the layers, (A^) — A^ -1 )) ^ 0. A similar phase 
called the a phase, where up-spin and down-spin Fermi 
surfaces are elongated perpendicular to each other, was 
reported^ 

We study the phase diagram as a function of bilayer 
coupling t± and chemical potential \i. The phase dia- 
gram shown in Fig. [2] is obtained for Fijlt = 0.8 and 
G2/F2 = 0.1. It shows that the parallel nematic phase, 
labeled by (1) and (2) in Fig. [2 is favored along the 
diagonal region (/i ~ t±) of the phase diagram. The 
hidden nematic phase (3) on the other hand, emerges 
as an intermediate phase at intermediate values for t± 
and is separated by a second order phase boundary from 
the neighboring isotropic regime (4). All the other phase 
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FIG. 2: (Color online) Phase diagram of the bilayer system as 
a function of inter- layer hopping t±/2t and chemical potential 
n/2t for F 2 /2t = 0.8, G 2 /F 2 = 0.1 and fa/fo = 0. The 
parallel nematic phase is labeled by // [regions (1) and (2)] 
and the hidden nematic phase by _L [region (3)]. The white 
areas denoted by (4) and (5) are regions of isotropic phase. 
All solid lines are first order boundaries, while the dashed line 
represents a second order line. The insets show representative 
Fermi surface topologies of the hybridized energy bands for 
different parts of the phase diagram. Note the different Fermi 
surface configurations for parallel order in regions (1) and (2). 
The phase diagram is symmetric in \i (due to particle-hole 
symmetry) and in t±. The DOS at various points (a)-(d) is 
shown in Fig. [3] and the conductivity and magnetization at 
points (h)-(j) are shown in Fig. [5] 

boundaries in Fig. [2] are of first order, and involve a sud- 
den change in the magnitude and/or orientation of the 
Fermi surface distortion dubbed nematicity. Note that, 
while the nematicity suddenly changes, there is no further 
symmetry breaking associated with the transition from 
one nematic phase to another. We call such a transition a 
"meta- nematic" transition in analogy to a meta-magnetic 
transition, where the magnetization jumps without any 
further symmetry breaking. The insets in Fig. [2] display 
the different Fermi surface topologies associated with the 
nematic phases in the bilayer system. 

To gain a better understanding of the different types 
of nematic order in the phase diagram, we study the be- 
havior of the DOS, since it was shown that nematic or- 
der develops in order to avoid a van Hove singularity^ 
In Fig. [3l we present the DOS for the points (a)-(d) as 
marked in Fig. [2l Originating from the underlying tight- 
binding dispersion, the bilayer DOS exhibits two peaks 
in the absence of nematic order, which are separated by 
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FIG. 3: DOS for the points (a) - (d) as marked in Fig. [2] (DOS in units of 2t/E). The parameters are fixed at F 2 /2t = 0.8, 
G2/F2 — 0.1 and t±/2t = 0.1. The order parameters have been determined self-consistently through Eq. JSJ). The formation 
of parallel nematic order splits the van Hove singularities into four peaks as shown in (b) and (c), while the hidden nematic 
order shifts the two peaks further apart from each other as shown in (d). 



2t± as shown in Fig. [3] (a). 



IV. EFFECT OF IN-PLANE MAGNETIC FIELD 



As we discussed above, nematic order prevents the 
Fermi level from lying at the van Hove singularity. How- 
ever, there is more than one channel to avoid a van Hove 
singularity in the bilayer system. Each of the original 
peaks can split into two new singularities, which leads 
to the parallel nematic phase and is similar to the single 
layer case. In total, this gives rise to four singularities 
as shown in Fig. [3] (b) and (c), where the separation be- 
tween two new peaks is given by 4 (1 + | A^ A ) | . In the 
hidden nematic phase, in contrast, the two original peaks 
shift further away from each other such that the mutual 
separation becomes greater than the bare hybridization 
2t± as shown in Fig. [3] (d). The peak separation then 

becomes 2 ^4(1 - |A( A )| 2 + (t ± ) 2 , which is always 

greater than 2t±. This channel, leading to the hidden 
nematic phase, is absent in a monolayer system. 



As we discussed above, the two consecutive meta- 
magnetic transitions and the large residual resistivity 
bounded by the metamagnetic transitions can be un- 
derstood within the nematic order proposal for a single 
layer square lattice. However, the recently discovered 
anisotropy between two longitudinal resistivities induced 
by tilting the magnetic field away from the c-axis cannot 
be explained within the nematic order scenario in a sin- 
gle layer, since the Zeeman coupling and the energetics of 
possible domains in a single layer system are independent 
of the magnetic field direction. 

However, note that Ru in Sr3Ru207 forms a bilayer 
layer square lattice, and the in-plane magnetic field has 
a dramatic effect in the bilayer lattice. For example, 
a field B x x leads to a k y mismatch between the up- 
per and lower layer, where the momentum difference 
p = 2ir/a (0 x /0o) y 5 while k x remains unchanged. The 
x-y symmetry breaking field lifts the degeneracy between 
the two nematic orientations, and thus the phase sepa- 
ration with domains is no longer energetically favorable. 
Hence the system recovers a pure nematic phase with 
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FIG. 4: (Color online) Phase diagram in the presence of a 
finite magnetic field (ty 0cc/0o = 0.2) as a function of the inter- 
layer hopping t± and chemical potential \i for F2/2t = 0.8 and 
G2/F2 = 0.1. No coupling of the magnetic field to the spins is 
assumed. All phase boundaries (solid lines) are of first order, 
as expected. The dashed lines indicate the phase boundaries 
in the absence of a field (cf. Fig. [2}. The magnetic field 
determines the direction of Fermi surface elongation. In the 
parallel phase (2), a rotation of the Fermi surfaces occurs 
across the white dashed line. The insets show the typical 
Fermi surface shapes encountered in various regions of the 
phase diagram. The DOS at the points (e)-(g) are plotted in 
Fig. E] 



an anisotropy in the longitudinal transport. Before we 
present signals of the anisotropy in various quantities, 
let us first study how the in-plane magnetic field affects 
the phase diagram of Fig. [2l The phase diagram in the 
presence of an in-plane magnetic field (n c/) x / $0 = 0.2) 
is shown in Fig. [U At this point we do not take into 
account the Zeeman coupling, but consider it in the sub- 
section below. The phase boundaries in the absence of an 
in-plane field (dashed lines) are also plotted, to make a 
comparison with the case without the in-plane magnetic 
field. 

Since the in-plane magnetic field breaks the x-y sym- 
metry, the hidden nematic phase denoted by (3) in Fig. [H 
is suppressed under the in-plane field. Naturally, the sec- 
ond order transition between the isotropic and hidden ne- 
matic phase in the absence of an in-plane magnetic field 
changes to a first order transition due to the presence 
of a symmetry breaking field. On the other hand, the 
regions with parallel nematic order are enhanced by the 
in-plane field, which is also expected. While the suppres- 



sion/enhancement of the hidden/parallel nematic phases 
under the in-plane magnetic field are rather robust fea- 
tures, the orientations of the Fermi surface elongations 
shown in the insets depend on the details of the bare 
band dispersion and the location of the van Hove singu- 
larities. In the Appendix, we present the DOS at the 
points (e)-(g) in Fig. [4] under the in-plane magnetic field 
to understand the relation between the nematic orienta- 
tions and the van Hove singularities. 

A. Conductivity and Magnetization 

In this section, we present conductivity, magnetiza- 
tion and susceptibility as a function of the magnetic field 
strength B x (B = B x x). We then propose a possible con- 
nection between the nematic phase and the phenomena 
reported in the bilayer ruthenate, Sr3Ru207. In addi- 
tion to the momentum mismatch due to the magnetic 
field discussed above, the Zeeman coupling is important 
to take into account, since it acts as a spin dependent 
chemical potential. 

The longitudinal conductivities are computed using the 
following standard Boltzmann equation 

^ = 2e 2 E ("S))^)^]' «» 

k,cr,i/=±l °^k,a 

where we set Tk, originating from impurity scattering, 
to be constant, while np is the Fermi-Dirac distribution 
function and stands for the Fermi velocity. The mag- 
netization and susceptibility are given by 

M= 9J ^E M<U) " (10) 
and 

The conductivity, magnetization, and susceptibility as 
a function of in-plane magnetic field are shown in Fig. [5] 
for F 2 /2t = 0.8 and G 2 /F 2 = 0.1. The points marked as 
(h), (i), and (j) in Fig. [2] show the values of the chemical 
potential and the bilayer coupling used for each panel 
from left to right, respectively. 

Let us first consider the case when the system is in 
the hidden nematic phase in the absence of a magnetic 
field [for instance, close to point (i) in Fig. [2]. As the 
magnetic field increases, the system will undergo several 
transitions. The isotropic-hidden nematic phase transi- 
tion at a weak field is hardly visible in the conductiv- 
ity and magnetization in Fig. [5] (i), although the emer- 
gence of the nematic order parameter is accompanied by 
an anomaly in the susceptibility as indicated by the left 
arrow in the susceptibility figure. At larger magnetic 
fields isotropic-parallel nematic transitions occur and are 
clearly signaled by the onset and offset of the conductiv- 
ity anisotropy and the jumps in the magnetization. 
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FIG. 5: (Color online) Longitudinal conductivities, magnetization and susceptibility as a function of magnetic field strength 
for i r 2/2t = 0.8, G2/F2 = 0.1 7/2t = 1.0. The three panels correspond to the points (h)-(j) as displayed in Fig. [2] The arrows 
in the susceptibility of the middle panel indicate the perpendicular-isotropic transitions. 



On the other hand, consider that the system is in the 
isotropic phase in the absence of a magnetic field, but 
close to the nematic instability [see e.g. points (h) and 
(j) in Fig. [2]. As the magnetic field is turned on, paral- 
lel nematic order develops in either the up or the down 
spin species, while the Fermi surface of the other spin 
species stays isotropic. Then the longitudinal conductiv- 
ity and the magnetization show clear indications of the 
isotropic-parallel nematic transition. While the preferred 
direction of the Fermi surface elongation in the presence 
of the symmetry breaking field depends on the details of 
the band structure and the chemical potential, a differ- 
ence between <j xx and <j yy is clearly visible in (h) and 
(j) in Fig. [5j Considering the in-plane field along the x 
direction, we find that the conductivity perpendicular to 
the field (<J yy ) is higher near half- filling as shown in Fig. [5] 
(h). Since the Hall conductivity <j xy is much smaller than 
the longitudinal components, the resistivity is almost in- 
versely proportional to the conductivity. Thus the resis- 
tivity parallel to the field direction is higher than that 
perpendicular to the field direction in panel (h), which is 
consistent with the observed anisotropy in S^R^Oy^ 

V. DISCUSSION AND SUMMARY 

Motivated by the recent experiment on the bilayer 
ruthenate^, we address the following question: can one 
explain the magnetoresistive anisotropy within the ne- 
matic order scenario which successfully describes both 
the metamagnetic transitions and the large residual 



resistivit y 12 ! 13 ? To understand the motivation of our 
study, one needs to recognize that the formation of ne- 
matic order in a single layer is insensitive to the direction 
of the magnetic field, and thus at first glance, the nematic 
theory cannot account for all the existing phenomena. 
However, here we show that the recently discovered mag- 
netoresistive anisotropy can be explained within the pic- 
ture of nematic order when the bilayer coupling is taken 
into account. It is essential to note that the in-plane 
magnetic field leads to a relative momentum mismatch 
between the layers through bilayer coupling. When the 
field is along one of the in-plane crystal axes, it breaks 
the x-y symmetry. Therefore, the degeneracy of two dif- 
ferent nematic orientations is no longer present, and the 
system recovers a pure nematic phase with anisotropic 
resistivities. 

To study the effects of an in-plane magnetic field on 
the nematic phases in a bilayer system, we first identify 
distinct nematic phases in the bilayer system. While the 
nematic phase always breaks the x-y symmetry in the sin- 
gle layer system, we find that there is another route to 
form a different nematic phase called the hidden nematic 
phase where the x-y symmetry is preserved. The hidden 
nematic and the isotropic phase are separated by a sec- 
ond order phase transition. While the x-y anisotropy is 
absent in both phases, the relative rotational symmetry 
between the layers is broken in the hidden nematic phase. 

The effect of the in-plane field is rather straightfor- 
ward when one recognizes the importance of the bilayer 
structure in Sr3Ru207. As discussed above, since the in- 
plane field is an x-y symmetry breaking field, there is no 



7 



spontaneously broken x-y symmetry in the presence of 
an in-plane magnetic field. It is clear that the two-fold 
degeneracy is no longer present, and domains cannot be 
formed under the same mechanism as described in Ref. 

one of the Fermi surface elongations is energetically 
preferred over the other for any small amount of in-plane 
field. Thus the system recovers the intrinsic anisotropy, 
unless it is in the hidden nematic phase. Based on the 
analysis of longitudinal conductivities and magnetic sus- 
ceptibilities, we propose that S^RuOr is close to the par- 
allel nematic instability in the absence of a magnetic field. 

One may question the validity of the mean field the- 
ory adopted in the current paper, since it is widely ac- 
cepted that Sr3Ru207 is a strongly correlated material 
with a putative quantum critical point. It is true that a 
mean field theory breaks down close to a quantum crit- 
ical point due to large fluctuations, and it is plausible 
that large fluctuations are important to determine an ef- 
fective Hamiltonian. However, one should note that we 
attempt to describe the ordered state and its first order 
transition to the isotropic phase. Since a mean field the- 
ory works quite well deep inside an ordered state due to 
negligible fluctuations, and the effects of fluctuation on 
a first order transition are not as important as for a sec- 
ond order transition, we argue that the mean field theory 
with nematic order qualitatively captures the phenomena 
discussed above. 

There are strong indications that a magnetic field 
tuned nematic phase bounded by isotropic regions exists 
in the bilayer ruthenate. In addition to the phenomena 
discussed above, a recent scanning tunneling microscopy 
experiment under a c axis field revealed the splitting 
of singularities in the local DOS across the metamag- 
netic transition^ However, further experiments to de- 
tect a direct Fermi surface anisotropy, such as scanning 
tunneling microscopy under in-plane magnetic fields, are 
highly desirable^ A microscopic mechanism for the for- 
mation of a nematic phase beyond the effective model 
Hamiltonian^, and the effects of disorder in relation to 
a putative quantum critical point are also important sub- 
jects for theoretical studies which we will address in the 
future. 
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VI. APPENDIX 

Here we present the DOS (Fig. [6j) in the presence 
of an in-plane magnetic field, to understand the corre- 
lation between the preferred direction of Fermi surface 
elongation and the locations of the van Hove singulari- 
ties. We set 7T(j) x /(j)o = 0.2, t±/2t = 0.075, F 2 /2t = 0.8, 
and G2/F2 = 0.1. The values of chemical potential used 
for the DOS correspond to those at the points (e)-(g) 
in Fig. [4j In the absence of an in-plane magnetic field, 
there are two peaks separated by the bilayer coupling, 
2t±. However, the in-plane magnetic field splits each of 
the peaks into two peaks such that there are four sin- 
gularities in the DOS. The top row shows the DOS of 
the underlying tight-binding bilayer model in the pres- 
ence of an in-plane magnetic field, where the magnetic 
field induces a splitting of the tight-binding singularities 
into asymmetric peaks. Here, we force the nematicity to 
be zero in order to see the effect of finite nematicity on 
the DOS, which is shown in the bottom row of Fig. [6l 
The three bottom panels represent three different ways of 
avoiding van Hove singularities in the presence of a finite 
magnetic field, which can be found by comparing the top 
and bottom panels for each case. In (e), a finite nematic 
order with two different orientations for each layer turns 
four peaks into two peaks such that the Fermi level is 
further away from the modified singularity. In (f) and 
(g), parallel nematic order not only splits the two peaks 
near the Fermi level further apart from each other, but 
also shifts the weight of the DOS between the singulari- 
ties. In all three cases, nematic order leads to a further 
separation between the Fermi level and the nearby van 
Hove singularity. This analysis helps us understand the 
orientations of the Fermi surface elongation for the par- 
ticular set of parameters used here, and one should bear 
in mind that the preferred direction of the Fermi surface 
distortion is sensitive to the parameters used in the ne- 
matic theory, since van Hove singularities depend on the 
details of the band structure. 
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